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We study numerically particle acceleration by the electric field induced near the horizon of a 
rotating supermassive (M ~ 10 9 — 1O 1O M0) black hole embedded in the magnetic field B. We find 
that acceleration of protons to energy E ~ 10 20 eV is possible only at extreme values of M and B. 
We also find that the acceleration is very inefficient and is accompanied by a broad band Me V- TeV 
radiation whose total power exceeds at least by a factor of 1000 the total power emitted in ultra-high 
energy cosmic rays (UHECR). This implies that if O(10) nearby quasar remnants were sources of 
proton events with energy E > 10 20 eV, then each quasar remnant would overshine e.g. the Crab 
nebula by more than two orders of magnitude in the TeV energy band. Recent TeV observations 
exclude this possibility. A model in which O(100) sources are situated at 100 — 1000 Mpc is not 
ruled out and can be experimentally tested by present TeV 7-ray telescopes. Such a model can 
explain the observed UHECR flux at moderate energies E ~ (4 — 5) x 10 19 eV. 



I. INTRODUCTION 

A conventional hypothesis of ultra-high energy cosmic 
ray (UHECR) acceleration in extragalactic astrophysical 
objects has two important consequences. First, it pre- 
dicts the Greisen-Zatsepin-Kuzmin (GZK) cutoff [lj in 
the spectrum of UHECR at energy of order 5 x 10 19 eV. 
Whether such a cutoff indeed exists in nature is cur- 
rently an open question Second, it implies that the 
observed highest-energy cosmic rays with E > 10 20 eV 
should come from within the GZK distance ~ 50 Mpc. 
Moreover, under plausible assumptions about extragalac- 
tic magnetic fields supported by recent simulations Q, 
the propagation of UHE protons over the GZK distance 
is rectilinear and the observed events should point back 
to their sources. While sub-GZK UHECR were found to 
correlate with BL Lacertae objects 0, @ , no significant 
correlations of cosmic rays with energies E 10 20 eV 
with nearby sources were found 0. 

In view of the latter problem, a question arises 
whether there exist UHECR accelerators which can pro- 
duce super-GZK protons and are quiet in the electro- 
magnetic (EM) channel. If such quiet accelerators ex- 
isted, they could explain the apparent absence of sources 
within ~ 50 Mpc in the direction of the highest-energy 
events. This idea was advocated, e.g., in Ref. where 
sources of UHE protons were associated with supermas- 
sive black holes in quiet galactic nuclei (the so-called 
"dead quasars"). However, it was pointed out that 
most of the energy available for particle acceleration in 
such an environment is spent for EM radiation by the 
accelerated particles. As a consequence, the flux of TeV 
7-rays produced by such an accelerator may be at a de- 
tectable level. 

Recent observations by HEGRA/AIROBICC O, MI- 
LAGRO [HI and TIBET P3 arrays improved substan- 
tially upper limits on the flux of 7-rays above 10 TeV 
from the point sources in the Northern hemisphere. This 
may exclude completely a possibility to explain the ob- 



served super-GZK cosmic rays by the acceleration near 
supermassive black holes. The purpose of this paper is 
to analyze this question quantitatively. To this end we 
study particle acceleration near the black hole horizon 
numerically. Following Refs. 0, we restrict ourselves 
to the case of protons. The case of heavy nuclei acceler- 
ation, propagation and detection is phcnomcnologically 
very different and requires separate consideration |13| . 
In particular, heavy nuclei can be easily desintegrated 
already at the stage of acceleration. 

We would like to stress that it is not our pupose to 
construct a realistic model of a compact UHE proton 
accelerator, but rather to find out whether quiet com- 
pact accelerators may exist, even if most favorable con- 
ditions for acceleration are provided. To this end we 
minimize the energy losses of accelerated particles by 
considering acceleration in ordered electromagnetic field 
and neglect all possible losses related to scattering of the 
accelerated particles off matter and radiation present in 
the acceleration site. However, we take into account, in 
a self-consistent way, the synchrotron/curvature radia- 
tion losses which are intrinsic to the acceleration pro- 
cess. Clearly, this approximation corresponds to most 
favorable conditions for particle acceleration. In realistic 
models the resulting particle energy will be smaller, while 
emitted EM power larger. Therefore, our results should 
be considered as a lower bound on the ratio of electro- 
magnetic to UHECR power of a cosmic ray accelerator 
based on a rotating supermassive black hole. 

We find that the flux produced by a nearby UHE 
proton accelerator of super-GZK cosmic rays in the en- 
ergy band E 7 > 10 TeV should be at least 100 - 1000 
times larger than that of the Crab nebula. The exis- 
tence of such sources is indeed excluded by recent obser- 
vations At the same time, the constraints on the 
sources of sub-GZK cosmic rays are weaker or absent, see 
Sect. VII for details. 

The paper is organized as follows. In Sect. II we 
describe our minimum loss model in more detail. In 
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Sect. Ill we present the analytical estimate and the 
numerical calculation of maximum particle energy. In 
Sect. IV the self-consistency constraints on the param- 
eters of this model are considered which arise from the 
requirement that there is no on-site e + e~ pair production 
caused by emitted radiation. In Sect. V the calculation 
of the EM luminosity of the accelerator is presented. In 
Sect. VI observational constraints are derived. Sect. VII 
contains the discussion of the results and concluding re- 
marks. 



II. THE MODEL 

The model we consider is based on a rotating super- 
massive black hole embedded in a uniform magnetic held. 
Because of the rotational drag of magnetic held lines, 
the electric field is generated and leads to acceleration 
of particles. In the absence of matter, the correspond- 
ing solution to the Einstein-Maxwell equations is known 
analytically at arbitrary inclination angle of black hole 
rotation axis with respect to the magnetic field 0, 0] . 
We assume low accretion rate and small matter and radi- 
ation density near the black hole, and neglect their back 
reaction on the EM and gravitational fields. We also ne- 
glect the effect of matter on propagation of accelerated 
protons. This corresponds to most favorable conditions 
for particle acceleration, and therefore leads to a maxi- 
mum proton energy and minimum EM radiation. 

The model has three parameters: the black hole mass 
M, the strength of the magnetic field B and the inclina- 
tion angle %■ We consider the maximally rotating black 
hole with rotation moment per unit mass a = M. This 
maximizes the strength of the rotation-induced electric 
field. For a given injection rate and geometry, the above 
parameters completely determine the trajectories of ac- 
celerated particles and, therefore, their final energies and 
the emitted radiation. We reconstruct particle trajecto- 
ries numerically keeping track of emitted radiation and 
taking into account its back reaction onto particle prop- 
agation. 

We assume that protons flow into the acceleration 
volume from the accretion disk which is situated at 
larger radii. We model this accretion by injecting non- 
relativistic particles uniformly over the sphere of the 
Schwarzschild radius i?s = 2GM, which is two times 
larger than the horizon of the maximally rotating black 
hole. We follow the trajectories of particles which prop- 
agate toward the horizon and then are expelled from the 
vicinity of the black hole with high energies. It turns 
out that such trajectories exist only if the inclination 
angle of magnetic field with respect to rotation axis is 
large enough, \ J> 10°. For smaller inclination angles 
all particles which propagate toward the horizon are fi- 
nally absorbed by the black hole. This means that the 
stationary regime in which particles accreted onto the 
black hole are subsequently accelerated and ejected with 
high energies exists only at x ^ 10° ■ In this regime the 



change of inclination angle and injection radius do not 
affect strongly the maximum energies of particles. The 
details of the particle acceleration and the description of 
the numerical procedures will be given elsewhere |13| . 



III. MAXIMUM PARTICLE ENERGIES IN THE 
STATIONARY REGIME 

In the abscence of matter and radiation backgrounds, 
particle energies are limited by the radiation loss intrinsic 
to the acceleration process. For a general electric and 
magnetic field configuration the energy loss in the ultra- 
relativistic limit is given by 0] 
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where m, e and v are particle mass, charge and veloc- 
ity, respectively. We use this equation in our numerical 
modelling to calculate the electromagnetic radiation pro- 
duced by the accelerated particles and to account for the 
back reaction of this radiation on particle trajectories. 

Before presenting the numerical results it is useful to 
summarize some simple qualitative estimates, see e.g. 
0, 0, 0] . Consider particle acceleration by a generic 
electromagnetic field obeying \E\ ~ \B\. If the energy 
losses can be neglected, energies of accelerated particles 
are estimated as 
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where we have assumed that the size of acceleration re- 
gion R is of order of the gravitational radius of the black 
hole R ~ 2GM . However, if the magnetic field strength 
is high, one cannot neglect the synchrotron/curvature en- 
ergy losses. 

If no special relative orientation of the three vectors E, 
B and v is assumed, Eq. JIJ becomes 
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which is the standard formula for the synchrotron en- 
ergy loss. Equating the rate of energy gain d£/dt = 
eE ~ eB to the rate of energy loss, one finds that in the 
synchrotron-loss-saturated regime the maximum energy 
is given by (see e.g. poj l 
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Here we have assumed that accelerated particles are pro- 
tons; for electrons the maximum energy is much smaller. 

The critical magnetic field strength at which the syn- 
chrotron energy loss becomes important can be found 
from the condition that the estimates @ and Q give 
the same result, 
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Here it is assumed that R ~ i?s- This critical field cor- 
responds to particle energy 
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which is a maximum energy attainable in the synchro- 
tron-loss-saturated regime for a given black hole mass. 

Acceleration is more efficient (the loss Q can be or- 
ders of magnitude smaller) in the special case when E, 
B and v are nearly aligned. These conditions may be 
approximately satisfied in some regions around the black 
hole. In this case particles follow closely the curved field 
lines and curvature radiaiton loss 
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becomes the main energy loss channel for high-energy 
particles. For an order-of magnitude estimate one can as- 
sume that the curvature scale of the magnetic field lines is 
of the order of the size of acceleration region, R ~ 2GM. 
This translates into the following maximum energy 
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which gives for protons 
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where we have set again R — Rg . The range of applica- 
bility of eqs.JSJ) and @ is given by the same condition 
B > B cr if 

In the numerical simulations, we injected protons uni- 
formly over the sphere surrounding the black hole. We 
disregarded trajectories which start at injection sphere 
and move outward. Among protons which approach the 
horizon and are then expelled to infinity, we selected 
those which have maximum final energy. For the black 
hole mass M = 10 10 Af Q , the dependence of this max- 
imum energy on the strength of the magnetic field is 
shown in Fig. 1 by red line (upper curve). For ener- 
gies of order 10 20 eV and higher, the numerically calcu- 
lated curve approaches the limit © which corresponds to 
the curvature-loss-saturated regime. The acceleration to 
these energies requires magnetic fields in excess of 10 4 G. 
The necessary magnetic field is even stronger for smaller 
black hole masses, cf. Eq.©. The maximum energies of 
protons do not depend strongly on the inclination angle 
in a wide range of x- 



IV. CONSTRAINTS FROM PAIR 
PRODUCTION 

There is an important self-consistency constraint which 
does not allow to increase B and M independently in or- 
der to reach higher energies. The reason is as follows. In 
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FIG. 1: The results of numerical calculation of maximum en- 
ergies of accelerated protons and accompanying 7-rays are 
shown by crosses (solid lines are fits to numerical data). The 
dashed line is the estimate JUJ of the proton energies in the 
curvature-loss dominated regime. The shaded region corre- 
sponds to the magnetic field strength exceeding the threshold 
of pair production by the curvature 7-rays. The black hole 
mass is M = 10 10 M P) . 



our model it was assumed that the acceleration proceeds 
in the vacuum. However, at a sufficiently strong mag- 
netic field, photons of curvature radiation may produce 
e + e~ pairs. Electrons and positrons will be accelerated 
in turn and produce more photons, which will again pro- 
duce e + e~ pairs, etc. The plasma created by this cascade 
will then neutralize the electric field and prevent further 
acceleration of particles. For consistency of the model we 
have to require that the cascade does not develop. 

Consider this process in more detail. The energy e 7 of 
the curvature photons in the regime when particle ener- 
gies are limited by curvature losses is estimated as 
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(cf. Eq. ||SJ|). Remarkably, the photon energy does not 
depend on the particle mass. This means that proton- 
originated and electron-originated photons have the same 
energy. Numerically, one has 
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If this energy is enough to produce more than one e + e~ 
pair within the acceleration site, the instability may de- 
velop. 

Therefore, for the stationary operating accelerator the 
mean free path d of a 7-ray in the background of a strong 
magnetic field (see |21|) has to be larger than the size of 
the acceleration region 
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This requirement leads to the condition 
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on the magnetic field in the vicinity of the horizon. 

In numerical calculation of proton trajectories, we kept 
track of the emitted photons. For given parameters of 
the accelerator, we determined the maximum photon en- 
ergy. The dependence of this energy on the magnetic field 
strength is shown in Fig. ^ by green line (lower curve). 
Substituting the calculated photon energy into Eq. Qlljl 
we can check whether the accelerator is in the station- 
ary regime. The shaded region in Fig. ^ corresponds to 
non-negligible pair production. The results of numerical 
calculation are in good agreement with Eq. I|12fl . 

From Fig. ^ we conclude that acceleration of protons 
to energies higher than 10 20 eV is marginally possible in 
a small region of parameter space (M, B). The magnetic 
field strength B should be close to the threshold of pair 
production. The black hole mass M should be larger 
than 1O 1O M0. Such black holes are rare. For example, 
in Ref it is found that supermassive black holes in 
AGNs range in 10 6 5 - 1O 1O ' 2 M with the mean mass 
being 10 8 ' 9 M Q . The list of nearby (within 40 Mpc) can- 
didates for quasar remnants |2^] does not contain black 
holes with masses above 5 x 1O 8 M0. 

Under reasonable assumptions about the black hole 
mass, the acceleration to energies above 10 20 eV is not 
possible in the stationary regime discussed above (no par- 
ticle production in acceleration volume). However, in 
spite of the fact that the accelerator can not operate per- 
manently, it is possible that UHECR are produced during 
"flares" , or short episodes of activity of the accelerator 
which arc interrupted by discharges. The natural dura- 
tion of one flare is about the time needed for the charge 
redistribution and neutralization of electric field in the 
acceleration volume to establish. This can be roughly 
estimated as the light-crossing time Tfl arc ~ Rs/c » 10 
hours for the 3 x 10 9 M© black hole. During such flares 
the electromagnetic luminosity of the accelerator would 
be much higher than the luminosity produced in the sta- 
tionary regime, because the electromagnetic flux is dom- 
inated by the radiation produced by e + e~ pairs whose 
number density is much higher than the density of initial 
protons. Since we are interested in the possibility of hav- 
ing a "quiet" UHE proton accelerator, we concentrate in 
the next Section on the case of stationary regime when 
the parameteres of the model are tuned to B rs 3 x 10 4 G, 
M w 10 10 M Q . 



V. ELECTROMAGNETIC LUMINOSITY OF 
THE ACCELERATOR 

It is clear from Fig. ^ that the acceleration of protons 
to energies above 10 20 eV proceeds in the curvature-loss- 
saturated regime. In this regime most of the work done 
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FIG. 2: Numerically calculated ratio of electromagnetic lumi- 
nosity of the accelerator to the luminosity emitted in particles 
with energies £ > 10 20 eV is shown by crosses, with the solid 
line being the fit to the numerical results. The shaded region 
corresponds to the magnetic field strength above the thresh- 
old of pair production by the curvature 7-rays. The black 
hole mass is M = 10 10 M©. 



by the accelerating electric field is spent on the emission 
of curvature readiation rather than on the increase of 
particle energy. The ratio of the dissipated energy to the 
final energy of proton is 
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Thus, the energy carried away by photons is at least hun- 
dred times higher than the energy carried by cosmic rays. 
Since only a small fraction of the accelerated protons 
reaches the UHECR energies £ > 10 20 eV, the ratio of 
the electromagnetic luminosity of the accelerator to its 
luminosity in UHECR with E > 10 20 eV is even higher. 

Numerically, we calculated this ratio as follows. We 
summed energies of those protons which were acceler- 
ated above 10 20 eV, and summed energy emitted in 
synchrotron/curvature radiaiton (including the radiation 
emitted by protons which did not acquire sufficient en- 
ergy while being expelled to infinity). We then took the 
ratio of the two sums. 

The results of numerical calculation of the ratio of elec- 
tromagnetic and UHECR luminosities in the stationary 
regime are shown in Fig.Elby crosses. Variations are due 
to fluctuations in the precise positions of the injection 
points 30]. The qualitative behavior of the numerical 
results is easy to understand. Close to the "threshold" 
£ = 10 20 eV, the accelerator emits finite power Lem but 
does not produce UHECR with £ > 10 20 eV. There- 
fore, the ratio Lem/Luhecr diverges at £ — ► 10 20 eV. 
If the magnetic field is large, the maximum energies of 
particles increase as well. However, the ratio 1Z also in- 
creases for each particle according to Eq. (|13fl . and so 
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does I/em/ Luhecr- The minimal value of I/em/^uhecr 
is reached at 8 « 1.5 x 10 20 eV. The numerically calcu- 
lated minimum of Lem/^uhecr is by a factor of 10 larger 
than the estimate (|13fl . 

When obtaining the results of Fig. [21 we have taken into 
account only the curvature radiation produced by pro- 
tons. For the magnetic field strength above the thresh- 
old of pair production ~ 3 x 10 G (and, correspond- 
ingly, £ max > 1.3 x 10 20 eV, see Fig. 0, our results give 
the lower bound on the electromagnetic luminosity of the 
compact accelerator. 

One can see from Fig. [2] that the electromagnetic lu- 
minosity of UHE proton accelerator based on a rotating 
supermassive black hole is 
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(14) 



Since the typical energy of photons of curvature radia- 
tion is about 10 TeV (see Eq. I|l()|l). the above relation 
implies that such a source of UHE protons would be much 
more powerful in the 10 TeV band than in the UHECR 
channel. 



VI. OBSERVATIONAL CONSTRAINTS 

The fact that production of UHE protons in the 
"quiet" accelerator is accompanied by emission of the 
TeV 7-ray flux enables us to put strong constraints on the 
possibility of existence of such accelerators in the nearby 
Universe. Following Ref. let us assume that there are 
about N ~ 10 nearby UHECR accelerators which are not 
further than Dqzk ~ 50 Mpc from the Earth. If these 
sources give the major contribution to the flux of cosmic 
rays above 10 20 eV 0, 
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We have seen in the previous section that the lower bound 
on the ratio of electromagnetic and UHECR luminosity 
of such a source is 1Z > 10 3 . This means that the elec- 
tromagnetic flux from the source should be 
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which implies the total luminosity larger than 

Lem ~ 10 44 erg/s (17) 

at the distance -Dqzk- This would be an extremely 
powerful source of TeV radiation. For comparison, the 
flux of the Crab nebula at energies above 15 TeV is 
■Fcrab ~ 10~ n erg cm" 2 s" 1 PH- Thus, the hypotheti- 
cal "quiet" cosmic ray sources which would explain the 



observed UHECR flux should be 100-1000 times brighter 
in the TeV band than the Crab nebula. 

The possibility of existence of persistent point sources 
of this type in the Northern hemisphere is excluded by the 
measurements of HEGRA AIROBICC Array [hJ and by 
MILAGRO experiment [UHl]. The upper limit on the 
energy flux from an undetected point source of ~ 15 TeV 
7-rays provided by HEGRA/AIROBICC group [Hj is at 
the level of Fhegra & (2-=-3) -Fbrab- Much tighter upper 
limit was published recently by the MILAGRO collabo- 
ration H3 Fmilagro (0.3 ™ 0.6) F Crab . 



VII. DISCUSSION 

The model of particle acceleration near the horizon of 
a supermassive black hole which we have discussed in this 
paper is based on a number of assumptions: maximum 
rotation moment of the black hole, low matter and radi- 
ation density in the acceleration volume, absence of back 
reaction of the accelerated particles and their radiation 
on the EM field, uniform magnetic field at large distance 
from the black hole. These assumptions have one com- 
mon feature: they facilitate acceleration to higher ener- 
gies and minimize losses (and, therefore, the produced ra- 
diation) . We have found that even in these idealized con- 
ditions the acceleration of protons to energy E = 10 20 eV 
requires extreme values of parameters, M ~ 10 10 Af Q and 
B ~ 3 x 10 4 G. Moreover, the acceleration is very in- 
efficient: the total power emitted in TeV gamma rays 
is 100 — 1000 times larger than in UHECR. In view of 
recent TeV observations, this rules out some UHECR 
models based on this mechanism of acceleration, e.g. the 
model of several nearby dormant galactic nuclei ("dead 
quasars" ) which was aimed to explain observed UHECR 
flux with energy E > 10 20 eV. 

In a more realistic case the above conditions may not 
be satisfied completely, and the acceleration of protons 
to energy E ~ 10 20 eV in the continuous regime may not 
be possible. The synchrotron losses due to the presence 
of random component -Brand of the magnetic field can be 
neglected if 
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where 1Z is given by l|13fl . This means that the presence 
of a tiny (1% level) random magnetic field will lead to the 
decrease of maximal energies of accelerated protons and 
increase of electromagnetic luminosity of the accelerator. 
Note that the synchrotron radiation is emitted in this 
case at energies 
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The power is still given by the Eq. (|14fl . 

Even if the strength of the random component of mag- 
netic field is as small as 10~ 5 B, for electrons, which 
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are inevitably present in the accelerator, the synchrotron 
losses will dominate over the curvature losses. The elec- 
tromagnetic power emitted by electrons will be in the 
100 MeV - 10 TeV energy band (see Eq. JTSJ). Assum- 
ing that the density of electrons is of the same order as 
the density of protons, one obtains the same estimate 
(|14fl for the 100 MeV luminosity of the accelerator. This 
means that such an accelerator is not only a powerful 
TeV source, but is also, in fact, an extremely powerful 
EGRET source. 

Even if the idealized conditions are realized in Nature, 
corresponding objects must be extremely rare. Thus, 
only a very small fraction of (active or quiet) galactic 
nuclei could be stationary sources of UHE protons with 
energies above 10 20 eV. 

If the parameters of the model are not precisely tuned 
to their optimal values, one expects the maximum en- 
ergies of accelerated protons to be somewhat below 
10 20 eV. It is therefore interesting to note that most of 
the correlations of UHECR with the BL Lacertae ob- 
jects come from the energy range (4 — 5) x 10 19 eV. The 
central engine of BL Lacs is thought to consist of a su- 
permassive black hole; it is possible that the acceleration 
mechanism considered above is operating in these objects 
|26j |. This mechanism may also operate in the centers of 
other galaxies which may possess (super)massive black 
holes, including our own Galaxy where it may be respon- 
sible for the production of cosmic rays of energies up to 
~ 10 18 eV [TSUI. 

The constraints from TeV observations are different 
in this case. First, cosmic rays of lower energies prop- 
agate over cosmological distances, so the UHECR flux 
is collected from a much larger volume, and the num- 
ber of sources may be larger. Correspondingly, the TeV 
luminosity of each source is smaller. Second, TeV radia- 
tion attenuates substantially over several hundred mega- 
parsecs. Third, at E < 10 20 eV the ratio Lem/Ljjhecr 
is smaller. Consider e.g. the case of O(100) sources lo- 
cated at z ~ 0.1 with typical maximal energy at accelera- 
tor £ ~ 5 x 10 19 eV. According to Fig.^ in this case the 
typical energy of produced 7-rays is ~ 4 TeV. The flux 
of 7-rays in this energy range is attenuated by a factor 
10 — 100, while according to Eq. ljT3)l . Lem/ Ljjhecr ~ 
50. Therefore, one may expect Feu ~ (0.01 — 0.1) fbrab 
for the TeV flux from each of these sources. This is within 
the range of accessibility of modern telescopes. For ex- 
ample, the TeV flux from the nearby (z = 0.047) BL Lac 
1ES 1959+650, which correlates with arrival directions 
of UHECR HI, is at the level of 0.06.Fcrab during 
quiet phase and rises up to 2.9Fo a b during flares. Sev- 
eral other BL Lacs which are confirmed TeV sources have 
fluxes w 0.03fbrab, see e.g. Ref. [2g|. 

This paper mainly concerns the stationary regime of 
acceleration when the acceleration volume is not polluted 
by the creation of e + e _ pairs. To ensure this condition 
we required that magnetic field does not exceed the crit- 
ical value (|12f> . If the magnetic field is larger, the ac- 



celeration by the mechanism considered here can only 
happen during flares which are interrupted by the cre- 
ation of e + e~ plasma and neutralization of the electric 
field as discussed in the end of Sect. II VI Although we do 
not have a quantitative model of a flare, some features 
of this regime and its consequences for UHECR produc- 
tion can be understood qualitatively. Since there is no 
constraint on the magnetic field in this regime, the maxi- 
mum energies of accelerated protons may exceed 10 20 eV. 
However, the efficiency of the acceleration during flares 
must be much lower than in the stationary case. First, as 
follows from Fig. |3 the ratio Lem/ -t'UHECR is larger at 
large B. Second, the dominant part of the EM radiation 
is produced by the created electrons and positrons whose 
number density exceeds by far the number density of pro- 
tons. Thus, one expects that the ratio £em/£uhecr for 
this sources is much larger than in Eq. i|14|) . 

An UHECR accelerators operating in the flaring 
regime would produce approximately constant UHECR 
flux at the Earth. The reason is the time delay of protons 
due to random deflections in the extragalactic magnetic 
fields. This delay is of the order of ~ 10 5 [a/l°] 2 yr for 
a source at 100 Mpc, where a is the deflection angle. 
Since the time scale of flares (light crossing time) is of 
the order of day(s), the variations of UHECR flux would 
average away. On the contrary, the TeV radiation from 
such a source would be highly variable, with powerful 
"TeV bursts" and average energy flux in TeV band ex- 
ceeding that in UHECR by a factor of 10 4 or higher. 
Note that there exist tight constraints on transient TeV 
sources: the energy flux of a TeV burst of duration 10 5 s 
has to be less than 10- 10 erg/cm 2 s - 10F Cra b dill. 
As in the case of stationary accelerator, this constraint 
excludes the possibility to explain observed UHECR flux 
by a few nearby proton accelerators operating in the flar- 
ing regime. The hypothesis of several hundred remote 
sources is not constrained by TeV observations. 

To summarize, the model of compact UHE proton ac- 
celerators which operate near the horizons of supermas- 
sive black holes in galactic nuclei can explain only the 
sub-GZK events. A large number (several hundred) of 
sources situated at cosmological distances is required. 
Production of UHECR in such sources may be associ- 
ated with the blazar-type activity, TeV 7-radiation be- 
ing an important signature of the model, testable by the 
existing 7-ray telescopes. 
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